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Abstract 

There is no indication so far on the spin of dark matter particles. We consider the possibility 
in this work that a spin-3/2 particle acts as dark matter. Employing the approach of effective 
field theory, we list all possible 4-fermion effective interactions between a pair of such fields 
and a pair of ordinary fermion fields. We investigate the implications of the proposal on the relic 
density, the antiproton to proton flux ratio in cosmic rays, and the elastic scattering off nuclei 
in direct detection. While the relic density and flux ratio are sensitive to all interactions albeit 
at different levels, the direct detection is only sensitive to a few of them. Using the observed 
data and experimental bounds, we set constraints on the relation of couplings and dark particle 
mass. In particular, we find that some mass ranges can already be excluded by jointly applying 
the observed relic density on the one side and the measured antiproton to proton flux ratio or 
the upper bounds from direct detection on the other. 
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1 Introduction 



There is compelling evidence from astronomical observations that the dominant component of 
matter in our universe is invisible, dubbed dark matter (DM). After years of efforts the evidence 
is still mainly restricted to the scope of gravitational effects. There are now many on-going 
or approved astronomical and laboratory projects that will hopefully reveal in the near future 
whether the dark matter is composed of particles or astrophysical objects like massive com- 
pact halo objects or it is not required at all because of modified Newton dynamics. For brief 
overviews on the current experimental status, see for instance the recent talks in Ref. flTJ . 

From the theoretical point of view we are not short of candidates if the dark matter turns 
out to be composed of particles. The most popular one is the lightest supersymmetric particle, 
perhaps the lightest neutralino, in supersymmetric models. There are also extensive discussions 
suggesting that the lightest Kaluza-Klein particle in extra dimension models @ or the lightest 
T-odd particle in little Higgs models with T parity [3] could act as DM, and so on. For detailed 
reviews, see for instance, Ref. EJ. Since the models on which the proposals mentioned are 
based are yet to be verified, it is important not to forget about other alternatives. In this context, 
effective field theory serves as a useful approach since one can focus on the interactions relevant 
to DM searches and parameterize unknown underlying dynamics in terms of effective couplings 
0, while leaving the dynamics to be identified in dedicated particle experiments like high 
energy colliders. 

The effective field theory approach has been widely employed to study the detection of a 
scalar, spin- 1/2 fermion and vector DM particle E El [HI [9l [IHl [III [121 [121 [131 [13 [III H71 [HI 
[T9l l20l . In this work we consider the possibility of a spin-3/2 particle as a DM candidate, 
and investigate its features in various detection experiments. Although we do not know yet a 
quantum field theory for such a particle of a higher spin that would be renormalizable in the 
conventional sense, this is by no means a reason to exclude its physical relevance: there are 
hadronic resonances of higher spin, e.g., A(1232) of spin 3/2, which play an important role 
in nuclear physics. Another example of spin-3/2 particles is the well-motivated hypothetical 
gravitino, which is a gauge particle associated with spontaneously broken localized supersym- 
metry. We are aware that even the effective interactions involving the A(1232) resonance are 
still controversial, see for instance j[2T| for a summary of the relevant issues. Our results in this 
work are nevertheless immune from such uncertainties since what we will need is not more than 
the Lorentz covariance from which the polarization sum of a spin-3/2 particle is constructed. 
We also know that the gravitino itself was previously suggested as a DM particle. It interacts 
with ordinary particles with an essentially gravitational strength, and the leading interactions 
at low energies are those involving a single gravitino field (and another superpartner); see Ref. 
|[22l for a practical introduction to the relevant issues. In contrast, what we will consider in the 
following are the effective interactions that contain a pair of spin-3/2 fields. Such interactions 
naturally preserve certain parity if the latter is required to stabilize the DM particle. 

The paper is organized as follows. In the next section we exhaust all possible 4-fermion 
effective interactions that involve a pair of spin-3/2 fields and a pair of ordinary spin- 1/2 fields. 
We then compute the cross sections for the annihilation and elastic scattering processes. These 
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results will be utilized in section 3 to compute the secondary particle fluxes in the cosmic rays, 
the relic density, and the effective cross sections between the DM particles and nucleons. Using 
the observational and laboratory data we set constraints on the couplings of the effective inter- 
actions as a function of the DM mass. It turns out some of the mass ranges can be excluded 
using the currently available data. We recapitulate our results in the last section. 

2 Effective interactions and cross sections 

Suppose the stability of the spin-3/2 particle is protected by certain parity, its leading effective 
interactions would involve a pair of it. Amongst the possible interactions that are relevant 
to their detection are those that couple to a pair of ordinary fermions. We therefore restrict 
ourselves in this work to the 4-fermion interactions. 

A free particle of spin 3/2 and mass M can be described by a field ^ that has the mixed 
transformation properties of a Dirac four-component field and of a four-component vector field 
||23l . Its equation of motion is 

(#-M)¥ M = 0, (1) 

augmented with the constraint, y^^n = 0- Multiplying eq (OQ) from the left by y" and applying 
the constraint one gets d^m^ = as a consequence. The wavefunction of such a particle satis- 
fies the same equations and thus has four independent components as desired, that correspond 
to the four spin states in its rest frame. As a matter of fact, the wavefunction for a particle 
with four-momentum p and helicity A, U^(p, A), can be constructed using the Clebsch-Gordan 
coefficients in terms of the ones for a spin-1/2 Dirac spinor u(p,s) of helicity s and a spin-1 po- 
larization e(p,K) of helicity K |[24|. Our following calculation will not depend on the explicit 
form of Up, but its polarization sum P V fi(p) = Ea U v {p,X)U^{p,X) which is evaluated to be 
(see for instance, Ref. j[22| upon correcting the sign of the M term) 

*Vv(p) = -V + M) (r MV (p) - \fT^ M )f) , (2) 

with T flv (p) = g^ v — p^ip v / p 2 and p 2 = M 2 . Note that the factor (p"+M) can be equally well 
put on the rightmost. The equation of motion and the constraint imply that 

fP^{p)=P^{p) 7 v = 0, 
P^ V (P)= P HV(P)P V = > 

W -M)P MV (p) = P^ v (p)(J>-M) = 0, (3) 

which may be employed to verify eq ©. The polarization sum Qjj, v (p) for the antiparticle and 
its relations can simply be obtained from the above by M — > —M. We will assume in this work 
that the spin-3/2 particle is of Dirac nature. For a Majorana-type particle, the amplitudes to be 
computed later are either multiplied by a factor of two or just vanish. 
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The leading 4-fermion interactions between a pair of ^ fields and a pair of ordinary 
fermion fields / are of dimension six. For simplicity we will not consider flavor-changing 
interactions. The independent Hermitian bilinears of / are 

(a) : //, fiy 5 f, Jytf, fy^y 5 f, fo^f. 

The bilinears of ^ are similar with the only difference in that they have two additional vector 
indices: 

Consider the self-contraction of a pair of indices in the list (?) . Since the interactions will be 
exploited in the cases where the spin-3/2 particles are on-shell, the constraint for the free field 
still applies. This means that it is not necessary to consider the contraction between the fields 
and the /matrices. For instance, after a little algebra, we find g ap *¥ aOpo^ p = — , which 

however was already covered in the list (i). The contraction is thus restricted to be between the 
two factors of the ^ field: 

(ii) : V«, y a iy 5 y a , V> a ypV a , ^ajpjs^, V a O pa V a , 

while further contraction produces nothing new. 

All possible interactions are exhausted by multiplying the terms in the list (a) with those 
in the list (?) and (ii) respectively and contracting remaining indices with the signature tensor 
g^ v or the totally antisymmetric tensor e^ v P a . Some of the terms so obtained can be removed 
as redundant. For instance, using £^ vpa o pa = —ilc^y*, (in our convention 75 = iy°y l y 2 y> 
and £ 0123 = +1), we have e a PP a ffV a o pa ^ p = -2ff¥ a y 5 W a , e aPpa fiy 5 f^ a o p(J ^p = 
—2fiy5fm a y$m a , and so on. The final list contains the following 14 operators: 

<...,4 = /c7 M v/(* M ^ v , * M 75* V ), fo^iysf^^v, ^rs^v); 
<..,8 = ff(^a^ a ,^aiy 5 ^ a ), fmf^a^^aiys^y, 
<..,12 = fr^af^^afys^JY^sf^af^^afYs^y, 

< 3 ,i4 = {fo^fJo^iy 5 f)^ a o^ a . (4) 
The effective interactions are summarized as 

^int = £L G i>/> (5) 

/ '=1 

where all couplings G\ are real and have the same dimensions as the Fermi constant Gf. 

To prepare for the numerical analysis in the next section, we display the cross sections 
for DM annihilation and elastic scattering off a nucleus. Since J£{ nt contains many possible 
interactions, it looks sensible to treat one interaction at a time. When *F in a specific model 
happens to interact in several ways with an ordinary fermion of a given flavor, one should 
sum coherently their contributions to a scattering amplitude. The spin-summed and -averaged 
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cross section in the center-of-mass frame for the annihilation process — > ff through the 
interaction , is 

f f (G f ) 2 s ls-4m 2 f 

of = * f ^^^JlMrn}/s,M>/s), (6) 

where s is the center-of-mass energy squared, rrif and M are respectively the masses of the final 
(/) and initial Q¥) particles, and N* = 1 (3) when / is a lepton (quark). The dimensionless 
functions A; for various operators are listed in the Appendix. 

In direct detection of dark matter one measures the recoil energy of nuclei that have been 
struck by a DM particle in the local halo. The event rate and energy deposited are determined 
by the cross section between the two. The calculation of the latter is a hard task, connecting mi- 
croscopic interactions of DM particles with quarks to effective interactions with nuclei through 
the intermediate chiral dynamics of nucleons, incurring uncertainties at each step; see the first 
article in Ref. for a review. The difficulty is alleviated to some extent by the fact that the 
collision is nonrelativistic. In this case the DM particles only feel the spin or mass of a nucleus 
0. In the remainder of this section we outline the procedure of this calculation relevant to our 
case and present the results for cross sections between *F and a nucleus. 

To begin with, one builds effective interactions between *F and nucleons from those for 
quarks shown in eq ©. One assumes that the Lorentz structures are not changed but the inter- 
action strengths get corrected by chiral dynamics. Since the nucleons in a nucleus can be treated 
nonrelativistic for the purpose here, the relevant nucleon bilinears (and thus quark bilinears) are 
restricted to the spin-independent (SI) and spin-dependent (SD) ones. In the approximation of 
zero momentum transfer, one then takes the diagonal matrix element of the nucleon bilinears 
for a static nucleus. The SI part essentially counts the numbers of the protons and neutrons in 
the nucleus while the SD part gets its main contribution from the unpaired nucleon spin. For 
a reasonable estimate of the nuclear matrix elements one has to appeal to nuclear models; in 
particular, when the momentum transfer must be taken into account, a form factor is necessary 
that will reduce the cross section with nuclei. 

Following the above procedure, the effective interactions in eq © are first classified into the 
SI and SD parts for nonrelativistic quarks: 

J2g = 5,(^)£G^,x = 5,6, (7) 

q 

JSf$ = B y (y)Y,GjV q ,y = 9,10, (8) 

q 

J§? Z SD = B*0P)£GfA*, z=l,...,4,ll,...,14, (9) 
q 

where S q = qq, V q = qfoq, A k q = q^y^q, and B X j(^¥) and B^^) are the *P bilinears without or 
with a free spatial index respectively. Although the nucleon matrix elements of the pseudoscalar 
quark bilinears do not vanish (see the second paper in Ref. [|25l ). the induced pseudoscalar 
nucleon bilinear has a nuclear matrix element that is suppressed by the velocity of nucleons. 



5 



Thus the operators ^7 g do not contribute at leading order to the scattering of W off a nucleus ./V 
as in the usual practice. The scattering amplitudes for the ^-N scattering are 



« 2m N f§B x (U), (10) 



*x,s 

« 2m N b y N B y (U), (11) 
4*1^(4)^(1/), (12) 

where m/v is the mass of the nucleus with atomic mass number A and charge Z, and (jfy p is 
the matrix element of the k-th component of the nuclear spin operator. The scalar SI effective 
coupling gets contributions from the protons and neutrons contained in N, 

f N = Zf p + (A-Z)f n , (13) 
f P (n) ~ L x ~mT Tq + 2J f ™ \ t G *-^> (14) 

q=u,a,s 1 q=c,b,t 1 

where m p ( n ,q) 1% tne proton (neutron, quark) mass. For light quarks the constants are 
related to the pion-nucleon sigma term ||25l , while for heavy quarks fp^ = 1 — £ /r£ 

q=u,d,s 

enter via the trace anomaly ll26l . The vector SI effective coupling is easiest to get since V q just 
counts the number of valence quarks when sandwiched between the nucleon states: 

b y N = Zbl + (A-Z)by v b p = 2G y i + G d y , W n = G u y + 2G d y . (15) 

Finally, the SD effective coupling is, 

S Z n = L<^, (16) 

K = ^[(S P Wq + (Sn}A q ], (17) 

where measures the fraction of the proton (neutron) spin carried by the quark q 11271 . 
and {S p ( n \) is the expectation value of the z-th component proton (neutron) spin operator in the 
nuclear state with the highest 11281 . 

The spin-summed and -averaged cross sections for *¥-N scattering at zero momentum trans- 
fer are, 

1 



where the amplitudes squared are evaluated in the standard manner 

£RJF = {2m N f N f- £ B x (U)Bl(U), (19) 

spins f spins 

= (^Nblf 1 - £ B y {U)Bl{U), (20) 

spins m spins 

£^fJp = (4m^)%(7 2V + l)i £ B k z (U)B?(U). (21) 

spins *F spins 
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The *F spin sums are computed using its polarization sum. In the nonrelativistic limit, this is 
facilitated by noting that only the spatial components are nonvanishing 

Hp) = M{% + l)(Sij + ^yirj). (22) 
The end results for the SI and SD cross sections due to various interactions are respectively 



J" 



K 



4 



^1,12,13 
°0 



K 



J N (J N +l)(g^ U,U ) 2 



20 20 80 

T'T'T 



(23) 
(24) 
(25) 



where /! = m^M/(m^+M) is the reduced mass for the ^-N system. That other operators do 
not contribute in the nonrelativistic limit can also be understood explicitly. Since Ufi(p,k) is 
built from u(p,s) and £»(/?, k), its static limit can be readily obtained. We have Uq(p,X) = 
either from p^U^,{p,X) = or by choosing physical polarizations with £o = 0. Independently 
of the Lorentz index in U^, the limits for a Dirac spinor bilinear also apply to U^, with the non- 
vanishing bilinears being restricted to Ua'fUp ~ O a Up, U a fy 5 Up, U a a' j Up « £ l i k U a y k, fUp, 
U a o Q, y s Up m iUaY^Up - Thus G\ A 6 7 8 vanish as they involve a pseudoscalar bilinear in 
or q, while 6\ disappears since it couples a nonvanishing q bilinear to ^o- Similarly, &f Q n 14 
do not contribute either at the leading order. 



3 Constraints from observations and experiments 
3.1 Relic density 

The dark matter produced in the early universe would either be depleted too much or over dense 
in the current epoch, depending on the interaction strengths in eq ©. The observed value for its 
relic density can therefore set constraints on the relevant parameters. To obtain the relic number 
density n^, one solves the Boltzmann equation 

^ + 3Hn ^ = _ (a | v |) _ (n ^2] > (26) 



where H = ySnp/SM^ is the Hubble expansion rate and is the value at thermal equilib- 
rium. Assuming the DM particles have negligible chemical potential we have n^> = nq, so that 
n T>M = 2n<p. (c|v|) is the thermally averaged annihilation cross section for a relative velocity v, 
and can be calculated in the reference frame where one of the *F particles is at rest ||29l . Using 
s = 2M 2 [l + (1 - v 2 ) -1 / 2 ] fa 4M 2 (1 + v 2 /4) for nonrelativistic DM particles, one expands eq 
© as, <j|v| = a + bv 2 + 0(v 4 ). Eq (l26l) is then solved numerically to yield [|30l . 

2 2xl.04xl0VGeV- 1 

1,2dM" ~ 1 ) 

M?\ a/ g* (xf) (a + 3b /x F ) 



7 



where g*(xp) is the number of relativistic degrees of freedom at the freeze-out temperature Tp, 
and Xf = M/Tf. The latter is solved self-consistently by 



Xf = In 



c(c + 2) 



45 gMM n (a + 6b/ x F ) 




(28) 



8 2^ y/jjJZPj 

where c is an order one parameter (we take c = 1/2), and g = 4 is the spin degrees of freedom 
of the *F particle.We employ the values of g* as a function of temperature T obtained in Ref. 
GO- 
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Figure 1: The current relic density (dashed curves) of the spin-3/2 DM, Qdm^ 2 > is predicted as 
a function of its mass for various interactions with couplings in scenarios I (left panels) and II 
(right panels) respectively. The horizontal solid band shows the observed range with the best-fit 
value, ^dm^ 2 = 0.1123 ±0.0035 021 . The number next to the legend indicates the operator 



We show in Fig. [j]the predicted current relic density for various interactions shown in eq © 
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as a function of the mass M. For simplicity, we consider two scenarios for the couplings G\ ifTTI . 
In scenario I, we assume a universal value for all interactions, g{ = 10~ 5 , 10~ 6 , 10~ 7 GeV 2 . 
For the operators involving a chirality-flip bilinear of ordinary fermions, i.e., those excluding 

tt?9 12, it is easy to imagine that they might be proportional to the mass of the involved fermion. 

We therefore study the scenario II in which g{ (1 GeV/m/) = 1(T 5 , 1(T 6 , 1(T 7 GeV 2 . For 
the purpose of comparison we include the results for ^9,. ..,12 in scenario II. The predicted den- 
sity decreases as the coupling G, (mass M) increases for a fixed mass (coupling). Also shown 
(horizontal band) is the range of the observed DM relic density, corresponding to the best-fit 
value, Qqm^ 2 = 0.1123 ±0.0035 11321 . Assuming that one of the interactions in eq © be re- 
sponsible for the observed relic density we plot in Fig. [2]fhe required couplings as a function of 
M. Since <y|v| increases with G ; and M, Gj has to decrease as M increases in order to match the 
observed relic density. 

The predicted relic density drops abruptly when a new annihilation channel is opened with 
increasing M. Similarly, for the relic density fixed to the observed value the required coupling 
Gi drops suddenly at each new threshold as M increases. This is especially obvious in scenario II 
at the tt threshold where the effect is significantly enhanced. Most curves fall into one of the two 
groups while the one corresponding to the operator £?f 2 stands alone. This arises from different 
behavior in their thermally averaged cross sections, (<r|v|) ~ a + b(v 2 ), where the coefficients a 
and b correspond to the s- and /?-wave annihilation respectively. While the operators ^1,3,5.7,10 
only give a b term, all others have both a and b terms. In addition, amongst the latter operators 
only &\2 has an a term that is proportional to m 2 , which explains its unique behavior in the 
figures. 




M (GeV) M (GeV) 



Figure 2: Assuming one of the interactions in eq. © produces the observed relic density, its 
coupling is shown as a function of the DM mass M for both scenarios I (left panel) and II (right). 

3.2 Direct detection 

The experimental results in direct detection of dark matter are conventionally presented in terms 
of cross sections on nucleons. We follow this practice to show our results on the cross sections 
oL / \ for nonrelativistic scattering of *F off a proton (neutron) due to various operators fff. 
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Figure 3: The spin-independent ^-proton cross sections (dashed curves) are plotted as a func- 
tion of M at different couplings for the scalar upper panels) and vector (Gg, lower pan- 
els) interactions. The left panels are for scenario I with g{ 9 = 10~ 5 , . . . , 10~ 8 GeV -2 , while 

thee right panels are for scenario II with Gg 9 x (1 GeV/m/) = 10~ 4 , . . . , 10~ 8 GeV -2 . The 
solid curves are the upper bounds from the experiments CDMS II (2010) P4ll . EDELWEISS-II 
(201 1) 03, XENON100 (201 1) (35|, XENON10 (2008) 03, and ZEPLIN-III (201 1) 01- 
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As we explained in the last section, out of many possible interactions there are only two types 
of them that contribute to the SI cross section and three types to the SD one. Gxy p i n \ are still 
given by eqs (I23I24I25I) for the SI and SD cases respectively, with the following substitutions: 
f 5 N f 5 p{n y b N b l( n y 8n ~+ L q G q z A p q M , J N -» 1/2, and jU Mm p{n) /(M + m p{n) ). For 
the chiral parameters related to the scalar SI matrix element we use the values in Ref. |[33l : 
= 0.020 ±0.004, fff = 0.026 ±0.005, /}f } = 0.1 18 ±0.062 for the proton, and fj£ = 

0.014 ±0.003,/^ = 0.036 ±0.008,/^ =0.1 18 ±0.062 for the neutron. For the nucleon spin 
fractions carried by quarks that are required in the SD matrix element, we assume the values in 



A£ = A^ 
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Figure 4: The spin-dependent ^-neutron cross sections (dashed curves) are plotted as a function 

f f 

of M for the axial vector G\ 2 and tensor ^13 interactions. The solid curves are the upper 



bounds from the experiments CDMS 11(2009) |39l, XENON 10 (2008) gOJ and ZEPLIN-III 
(201 1) 11381 . The notations for couplings are similar to Fig. [3j 

In Fig. [3] we display the SI ^-proton cross sections c^f p as a function of the mass M in 
both scenarios of couplings, while in Fig. |4] we plot the SD ^-neutron cross sections o\, n — 
o^ n and o^ n . In scenario II, both light and heavy quarks contribute equally to the scalar SI 
Oy p f n \ while for the vector SI ^ p i n \ b°th contributions are suppressed by either a light quark 
mass or a vanishingly small content of heavy quarks in the nucleon. In the same scenario, the 
contributions to the SD o^,\ from both light and heavy quarks are significantly suppressed 
by either light quark masses or tiny spin fractions of heavy quarks in the nucleon. In scenario I 
with universal couplings both SI and SD cross sections are dominated by light quarks. Since the 
current experimental bounds on the SD cross sections are several orders of magnitude weaker 
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than the SI ones, the upper bounds that one can get on the SD effective couplings are also much 
weaker. 

3.3 Indirect detection 

The DM particles in our Galaxy can annihilate through the interactions in eq © to produce 
leptons and quarks that fragment and interact further with the interstellar gas to create more 
secondaries, including gamma rays, neutrinos, positrons and antiprotons. By comparing the ob- 
served cosmic ray fluxes with known astrophysical sources, it is possible to infer the properties 
of DM particles in our galactic halo and constrain their annihilation rate. We have employed the 
public computer code GALPROP pT| to simulate the antiproton to proton flux ratio. The code 
solves numerically with appropriate boundary conditions the transport equation for the number 
density of cosmic particles that takes into account diffusion and convection effects amongst 
others. For our purpose here, the source term in the equation will contain a piece due to the ^iPV 
annihilation 

where the sum is over all channels of quark production, (dN/dE)q is the antiproton number 
per unit energy produced in the qq channel, and p (r) is the mass density distribution of the DM 
particles. We use the Monte-Carlo program PYTHIA H2l to simulate the (dN/dE)q spectrum. 
In our numerical analysis, we take the NFW profile ll43l : 

P(r) _ r e 
P© r 

where p is the DM density at the solar location, r the distance of the sun to the Galactic 
center, and R the scale radius. We adopt the following values for these parameters from Table 
3 in Ref.|44l : p = 0.389 GeV cm" 3 , r & = 8.28 kpc, and R = 20 kpc. The Galactic DM 
particles should follow the Maxwell-Boltzmann velocity distribution. We choose the velocity 
dispersion v = a/ (v 2 (r )) = a/3/2v c (> ) with v c (r ) = 243.75kms~ 1 being the local circular 
velocity EH, so that (v 2 ) = 2(v 2 (r )}. 

In the calculation of the p/p flux ratio with GALPROP, the diffusion region of cosmic rays 
is described by a thick disk of thickness 2L pa 8 kpc and radius Rpa20 kpc, with the thin galactic 
disk of thickness 2h pa 200 pc and radius R lying in the middle. The charged particles traversing 
the solar system are affected by the solar wind, which results in a shift in the spectrum observed 
at the Earth compared to the interstellar one ll45l l46l . We have scanned the solar modulation 
potential 4> from 300 to 1000 MV, and found that 4> = 330 MV yields the minimal % 2 for the 
background flux. 

Although the PAMELA data on the p/p flux ratio can be accounted for by GALPROP 
based on the conventional propagation model of cosmic rays, it cannot exclude a small portion 
of contribution from DM annihilations. This will set a stringent bound on the annihilation cross 
section (c|v|). By varying it within the acceptable deviation ranges of the PAMELA data and 



I + Tq/R 

l + r/R 



(30) 
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evaluating the % 2 value, we obtain the 3a upper bounds on the couplings g\ for a given value 
of the mass M. The results are shown in Fig. [5]for both scenarios I and H. The behavior of the 
curves is quite similar to that shown in Fig. [2l but the drop is less steep as M increases. 

We do not consider here the PAMELA positron fraction excess and related effective in- 
teractions for a few reasons. It has been proposed that the excess could originate from some 
astrophysical sources such as supernova remnants or nearby pulsars that were not accounted for 
earlier, see for instance Ref. ||48l for a status review. If the excess is due mainly to the dark 
matter annihilation, a strong tension arises between the excess and the relic density that is ten- 
tatively parameterized by a 'boost factor' as large as a few hundreds or even a thousand, whose 
origin however is unclear ||49ll . Thus the excess itself cannot yet result in useful constraints on 
interactions with leptons. And finally we want to work out combined constraints in the next 
section where only the interactions with quarks are relevant in direct detections. 




M (GeV) M (GeV) 

f 

Figure 5: The 3(7 upper bounds on various couplings G- as a function of M as imposed by the 
PAMELA p/p spectrum ll47l . in scenarios I (left panel) and II (right). 

3.4 Combined constraints 

In this section, we present the combined constraints from the observed relic density, direct and 
indirect detection data discussed in previous subsections. Since the direct detection is only 
sensitive to the operators ^1,5,9,12,13, we show our results for these types of interactions. In Fig 
[6l we show the combined constraints on their couplings g{ 5 9 12 13 as a function of mass M in 
both scenarios I (left panels) and II (right ones). For the SI types of interactions ^5 9 the direct 
detection generally imposes a stronger constraint than the indirect detection except for in 
scenario II. For the SD types of interactions the direct detection is advantageous in some cases 
while the indirect detection is better in others. 

When the upper bound curves from direct detection (SI and SD) and PAMELA p/p data are 
located below the relic-density allowed curves, they can provide more stringent constraints and 
be used to exclude some regions of parameters for various types of interactions. In the scenario 
I of universal couplings, the PAMELA data excludes respectively the operators ^12,13 in the 
mass ranges of (10, 17) GeV and (10,47) GeV for any couplings. While the SI direct detection 
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Figure 6: Combined constraints on the couplings G\ 5,9, 1243 are obtained as a function of M 
from the observed relic density, direct detection experiments of XENON 10 (SD M'- neutron 
scattering) and XENON100 (SI ^-proton scattering), and the observed PAMELA p/p flux 
ratio. Left panels for scenario I and right ones for scenario II. 
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completely excludes the operators ^59, the SD one dominates for the operator &\ and excludes 
it in the mass range (12, 37) GeV. The situation for scenario II is similar but less stronger. For 
instance, G\ is not sensitive to either the SD direct detection or the PAMELA data so that any 
mass would be allowed, while only a smaller mass range (10,49) GeV for the operator ffq is 
excluded by PAMELA. 

4 Conclusion 

We have considered the option that a spin 3/2 particle acts as dark matter, and investigated 
the constraints on it imposed by the current observations and experiments. We worked in the 
approach of effective field theory and wrote down all possible 4-fermion effective interactions 
that involve a pair of spin-3/2 DM fields and a pair of ordinary fermion fields. Assuming one 
interaction at a time is responsible for the dark matter, we studied its implications on the relic 
density, the antiproton to proton flux ratio p/p in cosmic rays, and the elastic dark matter scat- 
tering off nuclei in direct detection. While the relic density and the flux ratio p/p are virtually 
sensitive to all interactions at different levels and for different scenarios of couplings, only a few 
are relevant to the direct detection experiments. These observational and experimental results 
can be employed in a complementary manner. Using the observed relic density one can predict 
the relation between the effective coupling and the DM mass for a given interaction. When 
this relation curve lies above the upper bounds set by the direct or indirect detection, we can 
exclude some of the parameter regions for the DM particle. For example, the SI XENON 100 
data excludes the whole mass range that we studied for the interactions 9 when the couplings 
are flavor universal. Depending on the types of interactions and scenarios of couplings, the SD 
direct detection data and measured flux ratio p/p can also exclude portions of mass ranges. 
Further precise measurements will help narrow down the survival windows thus far. 
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Appendix Functions Ai(r,R) 

The functions Aj(r,R) with r = m 2 /s and R = M 2 /s appearing in the annihilation cross 
section are obtained upon finishing the phase space integration: 
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